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We examine the spectral properties of three-dimensional quantum billiards with a single pointlike 
scatterer inside. It is found that the spectrum shows chaotic (random-matrix-like) characteristics 
when the inverse of the formal strength v~^ is within a band whose width increases parabolically as 
a function of the energy. This implies that the spectrum becomes random-matrix-like at very high 
energy irrespective to the value of the formal strength. The predictions are confirmed by numerical 
experiments with a rectangular box. The findings for a pointlike scatterer are applied to the case 
for a small but finite-size impurity. We clarify the proper procedure for its zero-size limit which 
involves non-trivial divergence. The previously known results in one and two-dimensional quantum 
billiards with small impurities inside are also reviewed from the present perspective. 
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I. INTRODUCTION 



The quantum billiard with pointlike scatterers inside is 
a solvable model which still retains most of the interesting 
characteristics of non-integrable quantum physics. The 
problem is based on obvious physical motivations. The 
billiard is a natural idealization of the particle motion in 
bounded systems. The one-electron problem in quantum 
dots is a possible setting which may be used as a single- 
electron memory, a promising computational device in 
the future. It is now possible to actually construct such 
structures with extremely pure semiconductors thanks to 
the rapid progress in the mesoscopic technology. How- 
ever, real systems are not free from impurities which af- 
fect the electron motion inside. In the presence of a small 
amount of contamination, even a single-electron problem 
in bounded regions becomes unmanageable. The model- 
ing of the impurities with pointlike scatterers is expected 
to make the problem easy to handle without changing 
essential physics, at least at low energy. 

In spite of its seeming simplicity, the billiard problem 
with pointlike scatterers is known to possess several non- 
trivial properties. In two-dimensional billiards with a 
single pointlike scatterer, one observes phase reversion 
of wave function with the adiabatic rotation of the scat- 
terer around certain points inside the billiard [y. This 
can be regarded as the simplest manifestation of the ge- 
ometrical phase or Berry phase [^jj^. Moreover, the two- 
dimensional quantum billiards with pointlike scatterers 
possess the properties of ultra-violet divergences, scale 
anomaly and asymptotic freedom which are analogous to 
the ones found in quantum field theories J0|. Also, there 
is a problem of so-called wave chaos [g|J6|; Through its 
wave-like nature, the quantum particle can be diffracted 
by pointlike scatterers, which should have no effect on 



the classical motion of the particle [[7|-|9| . 

A fundamental problem for quantum billiards with 
pointlike scatterers is to understand global behavior of 
the energy spectrum in the parameter space of particle 
energy z and the strength of the scatterers v. In particu- 
lar, statistical properties of the spectrum are important 
because they reflect the degree of complexity (regularity 
or chaos) of underlying dynamics. For two dimension, 
the problem has already been examined in details |l^,p| . 
The chaotic spectrum (level statistics of random-matrix 
theory Ilj-H) appears along the "logarithmic strip" in 
the parameter space {z,v^^). More precisely, the effects 
of a pointlike scatterer with formal strength v are most 
strongly observed in the eigenstates with an eigenvalue z 
such that 
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where M is the mass of a particle moviiig in the billiard 
and A is an arbitrary mass scale. Eq. (|l|) indicates that 
the maximal physical coupling is attained at the value 
of formal coupling that varies with the logarithmic de- 
pendence of the particle energy. This energy-dependence 
is a manifestation of a phenomenon known as the scale 
anomaly, or the quantum mechanical breaking of scale 
invariance [EsllTq]: In two dimension, the physics is ex- 
pected to be energy- independent, since the kinetic term 
(Laplacian) and the zero-range interaction (a ^-potential) 
are scaled in the same manner under a transformation of 
length scale. However, the quantization breaks a scale in- 
variance, and as a result, the strong coupling region shifts 
with a logarithmic dependence of energy. The condition 
(|l|) also shows that, for any value of formal strength v, 
the system approaches to the empty billiard when the en- 
ergy increases. Thus the system possesses the property 



of the asymptotic freedom. 

Quantum-mechanical biUiard problems with pointlike 
scatterers inside can be defined for spacial dimension 
d < 3. Contrary to the two-dimensional case, spectral 
properties in three dimension have been scarcely studied 
so far. The main purpose of this paper is precisely to 
fill this void. The logarithmic energy-dependence of the 
strong coupling region observed in two dimension has its 
origin in the energy independence of the average level 
density of the system. Since the level density is propor- 
tional to the square root of the energy in three dimension, 
one expects substantially different spectral properties. In 
this paper, we find that this is indeed the case. It is 
shown that the value of formal strength which induces 
the maximal coupling is independent of the particle en- 
ergy, whereas the width of the strip on which the strong 
coupling is attained broadens with square-root depen- 
dence of energy. This means that, in three dimension, 
for any v (7^ 0), the system exhibits chaotic spectra at 
the high energy limit. 

Another objective of this paper is to relate the findings 
in the purely pointlike scatterers to the realistic situation 
of small but finite-size impurities. For the pointlike scat- 
terers, the condition for the strong coupling also depends 
on the mass scale A which is introduced in the process of 
regularization. This reflects the fact that formal strength 
V does not have a direct relation to the observables as it 
stands. In order to clarify the physical meaning of v, 
we begin by approximating a finite-range potential with 
a (^-potential within a truncated basis. The size of the 
truncation depends on the range of potential. We then 
obtain a relation between the formal and bare strengths, 
the latter of which corresponds to the strength of the 
(5-potential within the truncation. The relation enables 
us to apply the results for pointlike scatterers to finite- 
range cases. Moreover, it clarifies the proper procedure 
and physical meaning of the zero-size limit of the finite- 
size potential in an intuitive fashion. 

The paper is organized as follows. In Sec. II, we de- 
duce, from a general perspective without any assumption 
on the shape of billiards, the strong coupling condition 
in three-dimensional billiards under which the effect of 
a pointlike scatterer becomes substantial. In Sec. Ill, 
we consider the case for a small but finite-size scatterer. 
By examining a relation between the formal strength of 
the scatterer and the energy eigenstates of finite-size po- 
tential, we rewrite the condition for the strong coupling 
in terms of the observables. The previously known re- 
sults for one and two dimensions are reviewed from the 
present point of view. We clarify the proper procedure 
and meaning of the zero-size limit of finite-size potential 
in one, two and three dimensions. We test the predic- 
tions in Sec. II with the numerical calculations in Sec. 
IV. We look at the level statistics of rectangular box with 
a single pointlike scatterer inside. In particular, the case 
where the scatterer is located at the center of the box is 
examined in details. We summarize the present work in 
Sec.V. 



II. CONDITION FOR STRONG COUPLING IN 
TERMS OF FORMAL STRENGTH 

Consider a quantum particle of mass M moving in a 
three-dimensional billiard of volume V. The eigenvalues 
and eigenfunctions of this system are denoted by En and 

- -^Vn{x) = Enipnix) (n = 1, 2, 3, • • •)• (2) 

We impose the Dirichlet boundary condition to the wave 
functions (/?„ at the billiard surface. The average level 
density at energy z has square-root energy dependence; 
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Suppose that a single pointlike scatterer is placed at 
X = xo inside the billiard. Despite the simplicity of 
a contact interaction, the Schrodinger equation suffers 
from short-distance singularities at the location of the 
scatterer, which needs to be renormalized. This can be 
done in most mathematically satisfying fashion in the 
framework of the self-adjoint extension theory of a sym- 
metric operator in functional analysis. Details are given 
elsewhere (see Ref. ]ll|). We just present the relevant 
results. Starting with the formulation of Zorbas ||l^, we 
obtain the equation for the eigenvalues of the system, Zn 
[n = 1,2,3, • • •), as 



G{z) = V- 



where 
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In Eq. (H) , v is the formal strength of the pointlike scat- 
terer and A is an arbitrary mass scale that arises in the 
renormalization. The formal strength v does not have a 
direct relation to physical observables as it stands. Its 
relation to physical strength of the scatterer is discussed 
later in Sec. III. Here we just mention the following two 
points; (1) To ensure the self-adjointness of the Hamilto- 
nian for the system defined by Eq. (Q) , one has to take v 
to be independent of the energy, and (2) In the limit of 
"y — > 0, the system approaches the empty billiard. 

The second term of G{z) in Eq. (J^) is independent 
of the energy z. It plays an essential role in making 
the problem well-defined; The infinite series in Eq. (|5|) 
does not converge without the second term. For spa- 
cial dimension d > 4, the summation in Eq. (^ diverges. 
This reflects the fact that the billiard problem with point- 
like scatterers is not well-defined for d > 4 in quantum 
mechanics. Within any interval between two neighbor- 
ing unperturbed eigenvalues, G{z) is a monotonically de- 
creasing function that ranges over the whole real number. 
Therefore, Eq. (0) has a single solution on each interval. 



The eigenfunction corresponding to an eigenvalue z„ is 
written in terms of the Green's function of the empty 
biUiard as 
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This shows that if a perturbed eigenvalue Zn is close to an 
unperturbed one En (or En+i), then the corresponding 
eigenfunction V'n is not substantially different from (/3„ (or 
tfn+i)- Thus, the disturbance by a pointlike scatterer is 
restricted to eigenstates with an eigenvalue around which 
G{z) has an inflection point. This is because each inflec- 
tion point of G{z) is expected to appear, on average, 
around the midpoint on the interval between two neigh- 
boring unperturbed eigenvalues. Let (z, G(z)) be one of 
such inflection points of G{z); z ~ {Em + i?m+i)/2 for 
some TO. In this case, the contributions on G{z) from the 
terms with n c^ m cancel each other, and we can replace 
the summation in Eq. (0) by a principal integral with a 
high degree of accuracy; 
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where we have defined a continuous function g(z) which 
behaves like an interpolation of the inflection points of 
G{z). In Eq. (0), {(pnixo)"^) is the average value of 
fnixo)"^ among various n. For a generic position of 
the scatterer, one has {tpnix^Y) — ^/V. Notice that 
G{z) ~ g{z) is valid only around the inflection points of 
G{z). Using an elementary indefinite integral 
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for z > 0, we obtain 

G{z) 
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The first term in Eq. (ph, which depends on the energy 
z, disappears both at _E = and E = 00. As a result, the 
average value of G{z) at the inflection points is indepen- 
dent of the energy. Eq. (g) indicates that the maximal 
coupling of a pointlike scatterer is attained with the for- 
mal strength v which satisfies 



^3/2^1/2 
2^r 



(10) 




FIG. 1. Typical behavior of G[z) with mass scale 

A = 1 in Eq. (5) and its linearized version is shown 
as a function of z. The latter is obtained by ex- 
panding G(z) at its inflection point z on the interval 
Em < z < Em+i- The coordinate of the inflection point 
is given by {i,G{i)) ~ ((S„ -I- £m+i)/2, ~M^ /2-k). Strong 
coupling is attained when v~^ takes a value within the range 
of the linearized function. 

The "width" of the strong coupling region can be es- 
timated by considering a linearized eigenvalue equation. 
Expanding G{z) around z, we can rewrite Eq. (0) as 



G{z) + G'{z){z-z)~v-^ 



or 



G'{z){z~2) 
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(12) 



In order to ensure that the perturbed eigenvalue Zm is 
close to z, the range of w~^ -I- Af'^/^A^/^/27r has to be 
restricted to 
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(8) where the width A is defined by 
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This is nothing but the variance of the linearized G{z) 
on the interval between the two unperturbed eigenvalues 
just below and above z (see Fig. m. The width can 
be estimated by the average level density at the energy 
under consideration as follows; 



A(z~) 



E 



\ {i - En)' 



:Pav{z) ^ 



'2Pav{z) 1 



-(.«(.^o)^)E{(_i),^„(,)-rp 

= TT^ {iPn{xo)^) Pav{z) 
M3/2 ^ 



(15) 



We have implicitly assumed in Eq. (|lj) that the unper- 
turbed eigenvalues are distributed with a mean interval 
Paviz)~^ in the whole energy region. This assumption 



is quite satisfactory, since the denominator of G'(z) is of 
the order of (z — Enf' ^ indicating that the summation in 
Eq. ( p^ ) converges rapidly. 
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FIG. 2. Plot of v-^ = -M5/27r ± A(2)/2. The effects of 
a pointlike scatterer on the energy spectrum are expected to 
appear in the eigenstates in the region between both curves. 



We recognize from Eqs. (13) and (IT^) that the effects 
of a pointhke scatterer of formal strength v are substan- 
tial only in the eigenstates with eigenvalues z such that 
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The width A is proportional to the average level den- 
sity, and as a result, it broadens with square-root energy 
dependence (see Fig. ||). This can be understood from 
another perspective, by considering a scale transforma- 
tion of a heuristic Hamiltonian with a ^-potential; 
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(17) 



Although the Hamiltonian (fijy is not well-defined in case 
of spacial dimension d > 2, we proceed further for the 
moment and return to this point in Sec. III. Under a 
scale transformation x — > x/a, the Hamiltonian (O) is 
transformed to 
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Since the energy z scales as z — > a^z, the strength v 
which scales as v — > av must have square-root energy 
dependence, which explains Eq. (M). 

The findings in this section are summarized as follows; 

(1) For a three-dimensional billiard, the effect of a 
pointlike scatterer on spectral properties is maximal 
when the formal strength of the scatterer satisfies v~^ ~ 
— M^/^A^/^/27r, irrespective to the energy z. 

(2) The width A (or an allowable error in v~^ to look 
for the effect) increases with square-root energy depen- 
dence. 

From these two, we conclude: 

(3) For any value of formal strength (ti 7^ 0) , a pointlike 
scatterer tends to disturb a particle motion in billiards, 
as the particle energy increases; The physical strength 



increases proportional to the square-root of the energy. 
This makes a sharp contrast to the asymptotic freedom 
observed in two-dimensional billiards. 

Before closing this section, we give a few words on 
the shape of the billiard. Our implicit assumption for 
the shape is that the average level density of the empty 
billiard is dominated by the volume term, which has a 
square-root dependence on energy. The assumption is 
justified for a generic three-dimensional billiard which 
has the same order of length scale in each direction, ir- 
respective to a full detail of the shape of the billiard. If 
one length scale is substantially smaller than the other 
two, the surface term dominates the average level den- 
sity in the low energy region. As a result, the spectral 
property at low energy is expected to change with the 
logarithmic energy-dependence which is specific to two 
dimension [see Eq. (|l|)]. 



III. FORMAL, BARE AND EFFECTIVE 
STRENGTHS 

As stated in the previous section, v in Eq. (^ does 
not have a direct relation to physical observables as it 
stands. The main purpose of this section is to clarify 
the physical meaning of the formal strength v. For this 
end, we relate the formal strength to a strength defined 
through a more realistic potential with small but finite 
range. The relation makes it possible to apply the find- 
ings in the previous section to the finite-size impurities. 
The previously published results in two dimension pT| ] 
and the well-known elementary results in one dimension 
arc also reviewed from the present perspective. 

We first point out that the definition of the formal 
strength is not unique. Indeed, Eq. (^) is not a unique 
candidate for G(z); It can be defined by any convergent 
series for z ^ _E„ which has a form. 



G(z) = lim 
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(19) 



Here, /„ is an arbitrary quantity independent of the en- 
ergy z, whereas it may depends, in general, on En and 
ipnixg). The first term in the parenthesis on RHS in 
Eq. (Il^) docs not converge as iV ^ cx) in spacial dimen- 
sion d = 2, 3. This means that /„ should be taken as 
a counter-term which cancels the divergence of the first 
term. Once such a series {/«} is chosen, one can de- 
fine an equation, G{z) = 7, with an energy-independent 
constant 7. This gives a possible eigenvalue equation for 
the billiard with a pointlike scatterer of a certain fixed 
(energy-independent) coupling strength. It is obvious 
that, even with another choice of series, say {/«}, the 
same eigenvalue equation can be reproduced by shifting 



the value of 7 by 
is given by 
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(20) 



with an arbitrary real number A (7^ 0). This choice along 
with the definition v = l/j gives precisely the original 
eigenvalue problem Eqs. (^ and (||). Clearly, this v is 
a mathematical quantity whose physical interpretation is 
not immediately evident. 

To reveal the meaning of the formal strength v in 
Eq. (H), we begin by approximating low-energy spectra 
(eigenvalues and eigenfunctions) of a finite-range poten- 
tial by that of a zero-range interaction. Suppose that a 
small but finite-size scatterer of volume fl is located at 
X = xo inside a three-dimensional billiard of volume V. 
We describe the scatterer in terms of a potential which 
has a constant strength on a region fi; 



Uix) = 



Uq, xeft 
0, xev-n, 



(21) 



where the regions of the potential and the outer billiard 
are denoted by the same symbols as the volumes. We 
assume that the scatterer has the same order of size, say 
R, in each spacial direction, and also assume that the 
volume of the scatterer is substantially smaller than that 
of the outer billiard; fi ~ i?^ ^ F. In this case, the 



scatterer behaves as pointlike at low energy z <$^ E 



where i?7v(n) is estimated as 
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(22) 



Furthermore, the coupling of higher energy states than 
EN{n) to the low-energy states is weak, since wave func- 
tions with wavelength shorter than R oscillate within 
the scatterer. This means that the low-energy states 
{z <C i?Ar(si)) can be described by the Hamiltonian in 
terms of a (5-potential, Eq. (p^), with the coupling 
strength 
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(23) 



together with a basis truncated at E]^m\ . The truncation 
of basis is crucial for the present argument. In case of 
spacial dimension d > 2, a (S-potential is not well-defined 
in the full unperturbed basis. This is clear from the fact 
that the summation in Eq. (KJ) diverges in the limit of 
f] — > (hence as N{n) — > 00). The finiteness of the 
scatterer introduces an ultra-violet cut-off in a natural 
manner, and as a result, the low-energy spectra can be 
reproduced by the Hamiltonian ( |l7|) within a suitably 
truncated basis. 

In an analogy to the terminology of the field theories, 
we call the coupling v as the bare strength^ since it ap- 
pears as the coefficient of the (J-potential within a given 
truncation pi]. The bare strength v can be related to 



formal strength v as follows. Within the truncated ba- 
sis { (pn{x) ; n = 1,2, ...,7V(ri)}, the eigenvalues of the 
Hamiltonian ([l7|) are determined by 
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Inserting Eq. (p3) into Eq. (H) with Eq. (@), we obtain 
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Eq. (25) gives an exact relation between formal and 
bare strengths. In order to have further insight on Eq. 
(p5|), we take an average for ipnix^)^ among various n, 
\fn(xQ)^) ~ 1/V^, and replace the remaining summations 
on RHS by integrals. We then have 
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Using Eq. (|3|), along with elementary integrals 
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we can rewrite Eq. (Eq) as 

^3/2^1/2 ^3/2 
-~'^--'-^ ^rj^Pi'\^^E,,,,). (29) 

In Eq. (p7h, the first term on RHS is negligible in case 
oi z '^ E. Hence, at low energy z <C i?Ar(n), we have 
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The third term on RHS in Eq. (^) diverges as i?Ar(o) 
increases. This is exactly the same divergence which we 
observe in the summation in Eq. (E4h [or Eq. w^j] with 



(27) 
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opposite sign. This ensures the convergence oi G{z) in 
Eq. (0). Using Eq. (B2), we arrive at 



(31) 
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In order to reproduce a zero-range scatterer with a 
fixed formal strength v (7^ 0), the RHS of Eq. (31) needs 
to converge as n shrinks. This means that, for small 17, 
V should take a form. 



v{n) = 1/ - 
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(32) 



The first term in the denominator is a counter term that 
cancels the divergence of the third term on RHS in Eq. 

(1); 



c 
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[More precisely, C should be taken to cancel the diver- 
gence which appears in the summation in Eq. ( |24| ) .] The 
remnant quantity r(Q) in the denominator is a regular 
function which converges as O — > 0. In the zero-size 
limit, the finite-size scatterer shrinks into a pointlike one 
with formal strength 



r(0) 



M3/2A1/2 

2^ ' 



In terms of the potential height Uq, Eq. 
as 



Uo{n) = 1/ (-cf72/3 + r{n)n\ 

= 1/ (^-cn^/^ + r{o)n + o(f7)) 



(34) 



2^ is rewritten 



Since C is positive, we obtain 
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Eq. (^^ indicates that the potential has to be negative in 
the zero-size limit, irrespective to a form of r{fl). This is 
consistent with the fact that a pointlike scatterer with an 
arbitrary formal strength v (7^ 0) sustains a single eigen- 
state with an eigenvalue smaller than Ei. A seemingly 
plausible limit, fl — > along with keeping Uofl constant, 
is not allowable in the case of three dimension; Such a 
limit induces too strong a potential to define a quantum 
mechanical Hamiltonian for a pointlike scatterer. Notice 
that Eq. (36) does not exclude a possibility of strong 
repulsion t/o ^ on a small-size region ft ^ 0. Indeed, 
Eq. ( 13^ ) [or Eq. (|35|)] does not impose any restriction 
on V (or U) for any finite fl. As long as 51 is finite, one 
can reproduce even a strong repulsion by taking r{Q) as 
slightly larger than C/fl^/^. Such r{Q) is, in general, a 
very large positive quantity which diverges to -foo when 



fi shrinks into the zero-size together with positively fixed 

Combining the findings in the current and previous 
sections, we can deduce the condition for the strong cou- 
pling for a finite-size scatterer. Inserting Eq. (pG) or Eq. 
(|l|) into Eq. (|l|), we obtain 
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for z < EN(^n)- Eq. (|8|) indicates that v'^ ~ -C/Q^^^ 
is the condition for the strong coupling, and hence that 
the effects of a finite-size scatterer at low energy most 
strongly appear when it is weakly attractive, namely, 
when the bare strength v is slightly negative. In the 
zero-size limit, the condition (pq) is equivalent to 



|r(0)| 
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A(z) 
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Eq. ( |39| ) shows that it is the inverse of r(0) that repre- 
sents a direct measure of coupling strength of the zero- 
size limit of a scatterer. This naturally leads us to a 
definition of the effective strength of a pointlike scatterer 
by 



"eff 
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Using the effective strength Veff , we can rewrite Eq. (|3 
as 
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It can be observed from Eq. (|34| ) that, if the ori- 
gin of u~i-axis is shifted to the strong coupling value 
_7Vf3/2Ai/2/27r, the formal strength v is identical to Veff- 
We can also say that v~2 is a "distance" to the strong 

coupling value v^^ ~ —C/H,^'^, which is, in general, a 
large negative quantity for small il. Inserting Eqs. (pq), 
df) and (|3|) into Eq. (ph, we have 



1 



2i/2Af 



Uon 7r2f7i/3 



< 



Af3/2 
^372"^ 



(42) 



Eq. (B2) is the condition for the strong coupling in 
terms of the "observables" ; At low energy where a finite- 
size scatterer can be approximated by a pointlike one 
{z <^ l/Mil2/3)^ a particle of mass M moving in three- 
dimensional billiards is most strongly coupled to a finite- 
size (~ £2) scatterer of potential height Uq under the con- 
dition ([42|). [As seen from the arguments above, the co- 
efficients of fi~i/3 and i/i in Eq. (^2|) are not exact, but 
they are of the order of, or approximately, the values in 
Eq. (H).] 



The effective strength of a pointlike scatterer can be 
defined in two dimension in a similar manner. However, 
an energy-dependent correction is needed in this dimen- 
sion. One possible way to show this is to follow the ar- 
guments in the previous and present sections. Ref. |0] 
has taken this path. Instead, we here take an alternative 
manner which makes it easy to understand the origin of 
the energy dependence specific to two dimension. We be- 
gin by reexamining the condition for the strong coupling 
in three dimension, Eq. (jlq), in terms of the (5-potcntial 
with a truncated basis. We start by rewriting Eq. ( p^ ) 



as 



m 



< 



A{z) 



(43) 



where g{z) is defined in Eq. (0). Recall that it behaves 
like an interpolation of the inflection points of G(z) in 
Eq. (||). The energy dependence of g{z) is expected to 
be small. Indeed, we have g{z) ~ — M'^/^A^/^/27r from 
Eq. ( P), irrespective to the energy z. This indicates that 
Eq. (03|) is equivalent to the condition (p6|). Inserting 
Eqs. (0) and d29) into the condition (Esh , we obtain 



«-l - (ipn{xof)P 



PavjE) 
Z-E 



dE 



< 



A(z) 



(44) 



This is the condition for the eigenvalue equation Eq. (|2J) 
to have a solution z around some inflection point of LHS 
in Eq. (|^. Using Eqs.(||) and (|^, we have 



Af3/2 (3) 



< 



A(z) 



(45) 



Eq. (HI) is identical to Eq. (|^) for z < E^^Qy Notice 

that F{ (z,0) = 0, namely, the lower bound does not 
contribute on the principal integral. 

Let us now consider a two-dimensional analogue of the 
finite-range potential (|ll); It takes a constant value Uq 
on a finite-size region of area Q and zero everywhere else. 
The bare strength v is defined hy v = Uofl a.s in three 
dimension. Then, one can deduce an analogous formula 
to Eq. (^^ in two dimension. A crucial difference in two 
and three dimensions lies in the energy-dependence of the 
average level density. For the billiard with area S, it is 
given by pav = MS/2tt, according to the Weyl's formula. 
Since pav is independent of energy in two dimension, the 
analogue of Eq. ( [44| ) is estimated by 



Fi^\z,E) 



In 



dE 

z~E 
\z-E\ 



A 



instead of Eq. (|27|). Using (93„(a;o)^) 
position of the scatterer, we obtain 



(46) 
1/5 for a generic 



M 
2^ 



{Fi'\z,E^ 



(2), 



'Nin))-Fr{z,0) 



< 



(47) 



namely. 



M 


f'^T- 


^^E^iny 


-z\ 


<^ 


27r 


l A 


A 


/ 


^ 2 



(48) 



Here, the width A is estimated in a similar manner as in 
Eq. (|T|); 



A ~ TT {ipn{X0) )Pav ^ -;^ , 



(49) 



which is independent of the energy z. The condition (Hq) 
is identical to Eq. (51) in Ref. |11| , apart from a minor 
change in the definition of the width A in RHS. In two 
dimension, F^ (^jO) does not vanish and indeed has a 
logarithmic dependence on energy. This is the crucial 
difference from the two dimensional case. At low energy 



z<s:E 



N{n) 



l/Mfi, we have 



M 
2^ 



In- 



A 



ln(MAr2) 



< 



(50) 



Eq.(pO[) indicates that as ft shrinks, v should behave like 

v{n) = 1/ (— \n{MAn) + r{nyj , (51) 

where r(f2) is a regular function which converges in the 
zero-size limit, ft — > 0. The first term in the denomi- 
nator ensures that the logarithmic divergence disappears 
in Eq. @). Inserting Eq. (|l]) into Eq. @), we obtain 
in the zero-size limit, a two-dimensional analogue of Eq. 



r(0) 



M 




p. 




In 





'2tt 




A 



< 



A 



(52) 



This indicates that one can define the effective strength 
of the pointlike scatterer by 



.;e;/(z)^l/(r(0)-^ln^ 



Eq. 



now reads 



'eff{z)' 



-^<^ 



(53) 



(54) 



Eq. ( |54| ) with Eq. (|53| ) embodies the logarithmic strip of 
wave chaos that is the condition for the strong coupling 
in two dimension. By comparing this to Eq. (|l|), we 
obtain 



^(0). 



(55) 



The effective strength v^/f can be regarded as the direct 
measure of coupling strength of the scatterer, as in three 
dimension, and its inverse, w^f , is a 'distance' to a loga- 
rithmic curve of the strong coupling limit. The logarith- 
mic energy-dependence in Veff , exhibits a specific feature 



in two dimension. It conies from non- vanishing F-l (z, 0) 
which can be regarded as the origin of scale anomaly in 
a formalistic sense. Eq. ( |5l| ) shows 

Uo{Vl) = 1/ ( ^ \n{MKn) + r{n)n 

= 1/ (^HMAn)+r{0)n + o{n)] . (56) 



Hence, we obtain 






-0, 

— oo, 



as il 



(57) 



This is consistent with the fact that a pointlike scatterer 
supports a single eigenstate with an eigenvalue smaller 
than El, irrespective to the value of formal strength v 
(^ 0). The condition for the strong coupling in two di- 
mension is rewritten in a comparable form with exper- 
iment. Inserting v — Uofl as well as Eq.(|9|) into Eq. 
(BQ), we obtain 



1 



M 



UqVI 27r 



\n{zMn) 



< 



ttM 



(58) 



for z <^ l/Mft. An arbitrary mass scale A disappears 
from Eq. (|58|). Similarly to the three dimensional case, 
the effects of a finite-size scatterer at low energy z ^ 
1/AIfl appear most strongly when it is weakly attractive 

Our treatment is also applicable to one-dimensional 
case. We end this section by showing that all the stan- 
dard results in the elementary textbooks on quantum me- 
chanics for one dimensional ^-function is recovered in our 
formalism. In one dimension, one can define a (5-potential 
(of strength v) with the full unperturbed basis. The sum- 
mation on LHS in Eq. (p4) is convergent in the limit of 
N{Cl) — > cx), since the average level density is inversely 
proportional to square-root of energy; 



Pav{z) 



2i/27r ^ 



(59) 



for one-dimensional billiards with width L. The con- 
dition for the strong coupling is given by an equation 
formally identical to Eq. ([44|); 



./n z- E 



< 



A(z) 



where ((p„(fo)^ 

A(z) ~ 7r^((y9„(fo)^)p, 



1/ L and the width is given by 

7rMi/2 1 



21/2 



(60) 



(61) 



The principal integral in Eq. 
the use of 



can be estimated with 



f[^\z,E)^ 



[z- E)^ 



dE 



:ln 



^/^+^/E 



yi- 



(z > 0). (62) 



Since we have F^ (2^,0) = F{ (z,oo) — 0, we get 



< 



A(z) 



(63) 



Therefore, in one dimension, the strong coupling with a 
pointlike scatterer is attained when the bare strength v 
is large. The property is energy-independent (no scale 
anomaly). Since the width becomes narrow as the en- 
ergy increases, the effect of a pointlike scatterer with any 
(finite) bare strength disappears in the high energy limit. 
The bare strength v is identical to the effective strength 
v^ff in one dimension. They are related to the formal 
strength by 



''elf 



CX3 

n=l 



^(^o) 



En 



El + A2 



(64) 



In contrast to two and three dimensions, no divergent 
quantity appears in the definition of effective coupling. In 
analogy to the similar situation in quantum field theories, 
one might call this property of one-dimensional pointlike 
scatterer as super-renormalizability. A pointlike scatterer 
of bare strength v is obtained as the zero-size limit of 
a finite-range {Q) potential with height Uq = v/D, in a 
natural manner. In order to ensure u 7^ 0, C/q should 
behaves like 



Uo{n) = l/(r(f7)f7) 



(65) 



where r{Q) is regular in the zero-size limit. Since no 
singular term appears in v{fl)~^ at 17 — » limit, the usual 
zero-size limit, in which the product Uq^ is kept constant, 
is attained by keeping r(Jl) constant as fl varies. Thus, 
one obtains a pointlike object with the bare strength 



l/r(0). 



(66) 



We may conclude from the current perspective that it is 
an accidental fortune of super-renormalizability, that has 
enabled the simple formulation of the one-dimensional 
Dirac (S-function with a straightforward limiting proce- 
dure. 



IV. NUMERICAL EXAMPLE 

We have revealed, in Sec. II, the condition for the 
appearance of the effects of a pointlike scatterer in three- 
dimensional quantum billiards. It has been applied to 
the low-energy spectrum in case of a small but finite-size 
scatterer in Sec. III. In this section, the predictions are 
confirmed by examining statistical properties of quantum 



spectrum. We restrict ourselves to the exactly pointlike 
case. Even in this case, the numerical burden of handling 
very large number of basis states is quite heavy, and a 
smart trick is required to overcome it. 

We take a rectangular box as an outer billiard. We 
also assume the Dirichlet boundary condition such that 
wave functions vanish on the boundary. The mass scale 
is set to A = 1 in the following. Fixing the value of A 
makes all parameters dimensionless. The eigenvalues £"„ 
and eigenfunctions </?„(a?) in Eq. ^ are given by rear- 
ranging the triple- indexed eigenvalues and eigenfunctions 
in ascending order of energy; 



En 



^n^nyu, {x) = y 77 sni — sm — — sm — — , (67b) 

{n^ny,n^ = 1,2,3,- • •). 

The mass of a particle and the side lengths of the bil- 
liard are assumed to be M = 1/2 and {Ixjyjz) = 
(1.0471976,1.1862737,0.8049826), respectively. In this 
choice of the side lengths, the volume of the billiard is 
V ~ 1. We calculate G{z) on the interval between Em 
and Em+i by 



G{z) ~ iVnixo)^) 



1 



E 



m+2000 



Y^ ipn{xof ( — - 



z-E E^ + 1 

En 



Pav{E)dE 



En El + 1 



{Vn{xoY) (68) 

— ^ -^T^^ Pav{E)dE. 



E„ 



z-E E^ + 1 



When m < 2000, the first integral is discarded and the 
lower bound of the summation is replaced hy n = 1 in 
Eq. (|68|). The integral in Eq. ( |6^ ) is easily calculated by 
using Eq. (@). The approximation by Eq. ( pq ) serves to 
lessen numerical burden considerably, keeping a sufficient 
numerical accuracy. 

For a moment, we restrict ourselves to the case where 
the scatterer is placed at the center of the billiard. In 
this case, ((^„(a^o)^) = 'fnixoY = 8/V, which is eight 
times larger than the average value for generic cases. 
However, Eq. ( [iq ) is still valid, since only eighth of 
the whole unperturbed states, namely, that with even 
parity in each direction are affected by the scatterer, 
(rixny^nz = 1,3, 5,---). The solid curves in Fig. repre- 
sent v~^ = —M'^/'^/2tt± A(2;)/2. According to the con- 
dition (nq) (with A = 1), the effects of a pointlike scat- 
terer are expected to appear mainly in the eigenstates 
in the region between both curves. This is in fact the 
case as observed in Figs. || and ^, where the nearest- 
neighbor level spacing distribution P{S) is displayed for 
various non- negative values oiv^^ in two energy regions; 



2:100 ~ ^3100 in Fig. | and Z17000 
respectively. 

Fig.3 



^20000 in Fig. |, 



f^.-!) [)fi 



, (67a) nx) m 




p{s) n.fi 
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FIG. 3. The nearest-neighbor level spacing distribu- 

tion P{S) is shown for various values of v~^ in case of 
the scatterer being located at the center of the rectangular 
solid. The statistics are taken within the eigenvalues between 
zioo = 1307.95 and Z3100 = 12932.70. (The eigenvalues are 
numbered by taking into account only the eigenstates with 
even parity in each direction.) The solid (dotted) line is the 
Wigner (Poisson) distribution. 

We have numerically confirmed that the sign reversion 
of v~^ does not change the qualitative behavior of the 
distribution in both energy regions. Figs. and show 
that the case of v'^ = -M"^!"^ j^-n = -0.056269769 ~ 
is closest to the Wigner distribution (solid line). It is 
numerically observed that the second moment of P(S) is 
given by J^ P{S)S'^dS ~ 1.5 for v~^ — 0, irrespective 
to the energy. This indicates that P{S) is Wigner-like 
in the whole energy region for v~^ ~ — M^'^/27r. As 
v~^ increases, P{S) tends to approach the Poisson dis- 
tribution (dotted line). However, its rate depends on the 
energy. While P{S) becomes intermediate in shape be- 
tween the Poisson and Wigner distributions at v~^ ~ 10 
in Fig. ^, such distribution appears at v^^ ~ 30 in Fig. 
^. This can be easily understood from Fig. g; The value 
of -Mi/27r + A(z)/2 is 11.3 at zieoo = 8303.96, and 
25.7 at 218500 = 42508.80, respectively. These values can 



be considered as the upper bound of v~^ for inducing a 
Wigner-like shape in P{S) at each energy region. With 
-0"^ beyond the bound, the system is not substantially 
different from the empty billiard, and as a result, P{S) 
resembles to the Poisson distribution. 
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FIG. 4. Same as Fig. 3 except for the energy region 
between Znooo = 40184.77 and 220000 = 44767.02. 
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FIG. 5. The spectral rigidity A3(L) is shown for various 
values of v^^ in the energy region between zioo = 1307.95 
and 23100 = 12932.70. The scatterer is located at the center 
of the billiard. The solid (dotted) line is the prediction of 
random-matrix (Poisson) statistics. 

In Fig. ^ the spectral rigidity A3(L) is shown for var- 
ious values of v~^. The average is taken in the same 
energy region as in Fig. H. We can see the gradual shift 



to Poisson statistics (dotted line) as v~^ increases. Be- 
yond v~^ ~ 20, the value of A3(L) is close to the Poisson 
prediction, i/15. There still exists an appreciable dif- 
ference from random-matrix prediction (solid line) even 
for the strong-coupling limit {v~^ = —M'^^^/2tt ~ 0). A 
similar tendency has been reported in two-dimensional 
cases pffl. This can be understood from the fact that 
the range of the n-th perturbed eigenvalue is restricted 
to the region between n-th and {n + l)-th unperturbed 
ones in case of a single pointlike scatterer. As a result, 
the number of perturbed eigenstates on a certain energy 
interval does not differ largely from the number of unper- 
turbed ones in the same region. This restriction does not 
apply to the case of multiple number of pointlike scat- 
terers. We can therefore expect that the increase of the 
number of scatterers makes the energy spectrum more 
rigid. For two-dimensional rectangular billiard, a recent 
calculation corroborate this argument M. 



Fig.fi 



P{S) 



1 



ff-l 

I: 



FIG. 6. The nearest-neighbor level spacing distribution 
P{S) is shown for v~^ = in case of the scatterer being lo- 
cated at a generic position in the rectangular solid. The statis- 
tics are taken within the eigenvalues between 2100 = 415.81 
and 23100 = 3503.68. The solid (dotted) line is the Wigner 
(Poisson) distribution. 

Up to now, we have placed a pointlike scatterer at 
a specific position, namely, the center of the rectan- 
gular box. We next show the level statistics for the 
case of a generic location for the pointlike scatterer. In 
Fig. ra, we show the nearest-neighbor level spacing dis- 
tribution P{S) for a box with a scatterer located at 
xo = (0.5129731,0.5489658,0.3342914). The formal cou- 
pling is chosen to be v~^ = 0. Although a nearly maximal 
couphng is expected to be attained with this value of w~^, 
the level repulsion is rather weak and the observed P{S) 
is considerably different from the Wigner distribution. 
This can be understood by considering the state depen- 
dence of ^pn{xof. In case that the scatterer is placed at 
the center, the value of (/?„(i'o)^ is independent of the un- 
perturbed states; (/?„(xo)^ = 8/F for even parity states 
in each direction. This ensures a smooth change of the 
value of G{z) at the successive inflection points. For a 
generic position of the scatterer, however, the value of 
fn(xof' changes nearly at random as n varies, causing a 
considerable fluctuation of the inflection points of G{z). 
As a result, it frequently happens that successive unper- 



10 



turbed states are not substantially affected by the scat- 
terer even with the strong coupling value of the formal 
strength. It should be also noticed that, for the generic 
position of the scatterer, the width of strong coupling is 
substantially smaller than its average estimate given in 
Eq. (|lq). This can be understood as follows. Define the 
width for the n-th state by 



A„(z) = Tr^(pn{xo)'^Pav{z) 



(69) 



with z ~ z„ ~ En- Since i^„(xo)^ ranges from to 
8/V as n varies, the width A„(z) fluctuates between 
and 8A(z) for a generic xq. Since its average is given 
by A(z), successive appearance of the width should be 
smaller than A(z). This also explains why the coupling 
of the pointlike scatterer is rather weak for the generic 
case. [For the case that the scatterer is located at the 
center, we have A„(z) = A(z), irrespective to the un- 
perturbed states.] Clearly, a successive existence of the 
eigenstates unaffected by the scatterer is specific feature 
of a single-scatterer case. As the number of scatterers 
increases, such tendency disappears because only in very 
rare occurrence, none of the scatterers has a substantial 
influence on successive unperturbed eigenstates, as long 
as the coupling strength of each scatterer satisfies the 
condition (Jig). Again, for two-dimensional cases, Nu- 
merical results support this assertion B. 



V. CONCLUSION 

To conclude the paper, we summarize the findings in 
the previous sections. Eq. ( p^ ) in Sec. II is precisely the 
necessary condition for the appearance of wave chaos for 
three-dimensional pseudointegrable billiards with point- 
like scatterers. The condition is essentially different from 
that for two dimension. Whereas it is described by a 
logarithmically energy-dependent strip with an energy- 
independent width in two dimension, it is characterized 
by a parabola with a symmetric axis parallel to the en- 
ergy axis in three dimension. This implies that in three- 
dimensional billiards, the effect of the pointlike scatterer 
is stronger in the higher energy region. The numerical 
experiments using the rectangular box confirm the asser- 
tion that even a single pointlike scatterer brings about 
wave chaos under the predicted condition, although the 
precise amount of the effect depends on the location of 
the scatterer. 

Since the condition for wave chaos, Eq. (|l^) is de- 
scribed in terms of the formal strength v of the pointlike 
scatterer, it is not directly applicable to the case for re- 
alistic finite-size impurities. For this in mind, we have 
examined a relation between formal strength v of the 
pointlike scatterer and the bare strength v of the finite- 
size potential which is defined in a natural way as the 
product of height and volume of the constant potential 
on a finite-size region. The relation between v and v also 
makes it clear how one should take a zero-size limit to 



obtain a pointlike object with a given formal strength. It 
is shown that v~^ has a inverse cubic-root divergence in 
ri — > limit in three dimension. It is also shown that one 
can use a regular part oiv~^ [''i^) in Eq. (|3^)] as a direct 
measure of the coupling strength of a small scatterer. In 
other word, the inverse of the regular part corresponds 
to the effective strength of the scatterer. Since the coef- 
ficient of the singular part of v~^ is negative, wave chaos 
is expected to appear at low energy in case of weak at- 
traction. 

We have reviewed the previously known results in two 
dimension from the present perspective. Similarly to the 
three dimensional case, the inverse of the bare coupling 
v^^ has to contain a singularity as a function of the size 
of the scatterer, and the regular part of v~^ plays a cen- 
tral role in determining the effective coupling strength. 
There is a crucial difference, however. In two dimension, 
a logarithmically energy-dependent correction term is re- 
quired to define the effective strength. The existence of 
the energy-dependent term results in a peculiar feature 
for two-dimension, namely, the scale anomaly. Its origin 
is identified as the z-dependence of F^ '{z,E = 0) in Eq. 
(Eq) for two dimension. There is no corresponding term 
for three (and one) dimension. 

In a sense, the current work amounts to the search of 
a sensible zero-size limit of small obstacles in the quan- 
tum mechanics of general spacial dimension. Apart from 
the case of one dimension, where super-renormalizability 
guarantees the existence of trivial limit (S function) , one 
encounters a subtle balance of divergence and renormaliz- 
ability, which results in non-trivial properties of coupling 
strengths. We hope that we have persuaded the readers 
that the model of the billiards with pointlike scatterer is 
a valid, mathematically sound, and practically useful ide- 
alization of the quantum system with small impurities. 
We also hope that the predictions in this paper are to 
be checked through the experiments in the laboratories. 
In particular, Eqs. ( ^ ) and ( p8| ) for three and two di- 
mensions respectively, can be directly tested, since they 
are stated in an experimentally controllable form. Re- 
cent progress of microwave techniques with macroscopic 
devices ]T9[-pT[ offers a possible opportunity. 

Numerical computations have been performed on the 
HITAC MP5800 and S-3800 computers at the Computer 
Centre, the University of Tokyo. We thank Prof. Izumi 
Tsutsui for helpful discussions. 
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